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We report a highly unusual angular variation of the upper critical field {Hc2) in epitaxial superlat- 
tices CeCoIn5(n)/YbCoIn5(5), formed by alternating layers of n and 5 unit-cells thick heavy- fermion 
superconductor CeCoIns with strong Pauli effect and normal metal YbCoIns, respectively. For the 
n = 3 superlattice, Hc2{0) changes smoothly as a function of field angle 0. However, close to the 
superconducting transition temperature, Hc2{0) exhibits a cusp near parallel field {0 = 0°). This 
cusp behavior disappears for n = 4 and 5 superlattices. This sudden disappearance suggests the 
relative dominance of the orbital depairing effect in the n = 3 superlattice, which may be due to the 
suppression of the Pauli effect in a system with local inversion symmetry breaking. Taking together 
the temperature dependence of Hc2{0), our results suggest that some exotic superconducting states, 
including a helical superconducting state, might be realized at high magnetic field. 



PACS numbers: 74.25. Op, 74.25.0p, 81.15.Hi 

In the absence of time reversal symmetry or space in- 
version symmetry, the Fermi surface (FS) can often be 
split into portions with different spin structures. To 
stabilize superconductivity under such conditions where 
spin degeneracy is lifted, unconventional pairing of quasi- 
particles is needed, leading to exotic superconducting 
states very different from the conventional BCS pairing 
state of (k -k 1). Considering the situation of the 
broken time reversal symmetry alone, Fulde and Ferrell 
[l] , and Larkin and Ovchinnikov [2] proposed the pairing 
state of (k -k+q l) on Zeeman-split FS. This so-called 
FFLO pairing state leads to the modulation of the super- 
conducting order parameter in real space with the mod- 
ulation wavelength of the order of l/|q|. On the other 
hand, in the lack of space inversion symmetry, a Rashba- 
type spin-orbit coupling splits the FS into branches with 
spins of opposite rotation sense [3 . When the magnetic 
field is applied to such a system, pairing state with a 
finite center-of-mass momentum can also be realized, re- 
sulting in a helical superconducting state analogous to 
the FFLO phase. 

However, such exotic superconducting states have been 
poorly explored because of the lack of suitable materials. 
Recent advancement in heavy fermion thin film fabrica- 
tion technology [4, 5 has enabled the preparation of su- 
perlattices formed by alternate stacking of c-axis oriented 
CeCoIns and YbCoIns with atomic layer thicknesses. 
The large Fermi velocity mismatch across the interface 
between CeCoIns and YbCoIns significantly reduces the 
transmission probability of quasiparticles, thereby en- 
suring quasi-two-dimensional superconductivity confined 



within CeCoIus layers [G', T . This provides a unique op- 
portunity to explore the physics discussed above. This is 
because bulk CeCoIns with strong Pauli effect has been 
reported to host the FFLO phase at low temperatures 
and high magnetic field jSHTT]. In the superlattice, the 
electronic structure becomes two-dimensional, which is 
expected to stabilize the FFLO phase [8]. Moreover, the 
three-dimensional magnetic order p!2Hl5] . which is re- 
sponsible for the perplexing situation in the bulk sample, 
is expected to be strongly suppressed [4 due to negligibly 
small RKKY interaction between the adjacent CeCoIns 
block layers through the YbCoIus spacer. Furthermore, 
the importance of local inversion symmetry at the inter- 
face between CeCoIus and YbCoIns in the superlattice 
has recently been suggested to play a decisive role on su- 
perconducting properties, in particular when the thick- 
ness of CeCoIus is only a few unit-cells thick ^16j . 

In this Letter, we report the precise angular depen- 
dence of Hc2 of the CeCoIn5(n)/YbCoIn5(5) superlat- 
tices with n = 3, 4 and 5. We find that, for the n = 3 su- 
perlattice, Hc2{0) exhibits a highly unusual temperature 
evolution, which can not be explained for conventional 
thin films or superlattices. Using Ginzburg-Landau (GL) 
theory, we show that this is consistent with the scenario 
that the ultra-thin, two-dimensional CeCoIns enters an 
inhomogeneous phase at low temperatures, and the anal- 
ysis of the n dependence allows us to discover the impor- 
tance of space inversion symmetry. 

The CeCoIn5(n)/YbCoIn5(5) superlattices used for 
this work were grown using molecular beam epitaxial 
technique, where n layers of CeCoIus and 5 layers of 
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FIG. 1: (Color online) (a) Temperature dependence of re- 
sistivity in n = 3 and n = 4 superlattices. Tc, defined as the 
temperature where the resistance is 50% of the normal state 
resistance, is 0.99 K and 1.04 K for n = 3 and n = 4, respec- 
tively. The inset presents the schematic of the superlattice. 
(b) Field dependence of resistivity in the n = 3 superlattice 
at 0.25 K at selected angles. Hc2 is similarly defined as the 
50% point of the normal state resistance. The inset shows 
the arrangement of the experimental configuration and the 
definition of 0. 



YbCoIns were stacked alternately, typically repeated for 
30 - 60 times (Fig. [l^, inset). The growth details and the 
characterization of the superlattices are described else- 
where [5 . Resistivity measurements were performed us- 
ing a standard four-contact method, with an electrical 
current of 3.16 jaA flowing along the 6-axis of CeCoIus. 
The superlattice was mounted on the rotator of a dilu- 
tion fridge, allowing the magnetic field direction to vary 
in the ac plane, with ^ = 0° corresponds to || a (Fig. 
[T]3, inset). In addition to the n = 3 superlattice where 
Tc = 0.99 K, we also measured n = 4 and n = 5 su- 
perlattices (Tc = 1.04 K and 1.24 K, respectively) for 
comparison. 

Fig. [l^ shows the temperature dependence of resis- 
tivity for n = 3 and n = 4 superlattices, where zero 
resistivity is observed. Since (,c{0)^ the zero temperature 
coherence length perpendicular to the conducting plane, 
is estimated to be about 3 unit-cells thick [8 , orbital de- 
pairing effect would ordinarily be weakened in these films. 
In this regard, it would also be interesting to study the 
n = 2 superlattice. However, the n = 2 system does 
not exhibit a full superconducting transition. When the 
magnetic field is applied, clear transition to the normal 
state is recorded, enabling an accurate determination of 
Hc2 as a function of angle (Fig. ^p)- 

The angular dependence of Hc2 in an anisotropic bulk 
superconductor is smooth for all 6>, and can be described 
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FIG. 2: (Color online) (a) Angular dependence of Hc2 for 
the n = 3 superlattice. (b) Hc2{0) of the n = 3 superlattice 
at 0.16 K and 0.8 K, normalised at ifc2(0°) and ifc2(90°), 
clearly exhibit different curvatures (c) For reference, Hc2(0) 
of the n = 4 superlattice at 0.25 K and 0.9 K, normalised at 
ifc2(0°) and ifc2(90°), are also plotted 
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However, for a thin film with thickness d smaller than 
^c, Hc2 obeys the following equation first derived by Tin- 
kham 18 : 
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Therefore, in the thin film limit Hc2 is non-differentiable 
at 0° and follows a cusp-like dependence at small 0. 

We present the Hc2{0) data of the n = 3 superlattice 
in Fig. [2k. At 0.8 K, Hc2{0) shows a kink at 0°, and the 
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FIG. 3: (Color online) The Hc2{0) of the n = 3 super- 
lattice, replotted in an appropriate dimensionless form (see 
text). The solid lines are the fits to the data using the model 
described in the text. 
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slope of Hc2{0) increases monotonically as |^| decreases. 
This is a characteristic behaviour predicted by Eq. ([2|. 
On the contrary, at 0.16 K, Hc2{0) is a lot smoother at 
0°, which is closer to the behaviour described by Eq. Q. 
For intermediate temperatures, the curvature of Hc2{0) 
appears to fall in between the two limits governed by 
Eqs. ([T]) and ([2|, respectively. The distinct functional 
form of Hc2{0) can be readily compared when we plot 
the datasets at 0.8 K and 0.16 K on the same axes, with 
normalised Hc2 at 0° and 90° (Fig. [^Jd). Note, however, 
that a similar exercise on the Hc2(0) data of the n = 4 
superlattice does not reveal such a drastic difference (Fig. 

From the theoretical point of view, the existence of 
the FFLO phase can be discussed in the framework of 
GL theory by including the 'gradient' term /c(|na^^|^ + 
|n^^|^) in the GL functional. In the FFLO state, k 
becomes negative. Therefore the inclusion of the term 
lowers the total energy and guarantees the stability of the 
inhomogeneous phase. In the homogeneous BCS state, 
k takes a positive value and it can be shown using the 
modified GL theory [Il9l |20] that for > the slope of 
the critical field assumes the following form near ^ = 0°: 
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where Hp{T) is the line of the second order transition 
taking into account only the Pauli paramagnetic effect, 
ao is a numerical constant, and ga is related to g-factor 

via g = \ g'^cos^O + g^sin^^. Therefore, the rounding of 



FIG. 4: (Color online) (a) The contribution of the linear and 
quadratic components, characterized by a and /S respectively, 
as a function of temperature for the n = 3 superlattice (b) 
The ratio oia//3 for the n = 3 superlattice as a function of the 
reduced temperature. The data for n = 4 and n = 5 superlat- 
tices are also shown for comparison. Inset: The temperature 
dependence of \dHc2/dO\Q^Qo extracted experimentally (see 
text). 



the cusp at 6> = 0°, which is manifested by the diminish- 
ing of \dHc2/d6\Q^Qo^ can come from two sources - (i) a 
decreasing k which accompanies the entrance to an in- 
homogeneous state, or (ii) the suppression of Hp which 
signals the enhancement of the Pauli effect relative to the 
orbital effect. 

We now attempt to quantify the relative merits of Eqs. 
([T]) and Q in describing our Hc2 data. We replot the 
data in Fig. |2] in the 'dimensionless' form, as shown in 
Fig. [3j Using this representation, the graph would be lin- 
ear for a pure Tinkham behaviour (Eq. [2| and quadratic 
for the case of a simple, anisotropic bulk superconductor 
(Eq. [T]). However, it is apparent that the overall angular 
variation of Hc2 follows an intermediate behaviour. To 
move forward, we propose a description of our data using 
the model below: 
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With this model, which essentially quantifies the rela- 
tive contributions of Eqs. ([T]) and ([2| through the coefii- 
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FIG. 5: (Color online) Schematic diagrams showing the im- 
portance of local inversion symmetry breaking with reducing 
n. The CeCoIns layers, shown in red, which possess (break) 
the local inversion symmetry are denoted by the open circles 
(crosses). With broken inversion symmetry, the FS is split 
into branches with spins of opposite rotation sense, and the 
degree of splitting is characterized by an. Note that some 
inner layers, e.g. the layers indicated by the closed circles 
for n = 5, also experience the effect of broken inversion sym- 
metry, albeit with a much weaker effect and hence a much 
reduced an. 

cients a and /3, an excellent description of the Hc2{0) is 
achieved, as evidenced by the solid lines in Fig. [3] 

The coefficients a and (3 extracted using Eq. Q are 
plotted in Fig. [4^ for the n = 3 superlattice. We find 
that |a| decreases as the temperature is lowered, whereas 
1/3 1 shows an opposite trend. We note that |a| and 
must sum to 1, as required by the construction of Eq. 
Q. This analysis confirms our visual inspection earlier 
that Hc2{0) follows the Tinkham-like behaviour closer at 
high temperatures. To gauge the relative contribution 
of a and /3, we calculate the ratio As illustrated 

in Fig. [4]3, for the n = 3 superlattice, this ratio de- 
creases drastically when temperature decreases. Hence, 
at low temperatures, the squared term described by P 
overwhelms the linear term associated with a. The con- 
trasting behaviour of the n = 4 and n = 5 superlattices 
is consistent with a much weaker temperature variation 
of \a/p\. 

The magnitude of a is directly proportional to 
\dHc2/dO\Q^Qo , since in the presence of the cusp both \a\ 
and \dHc2/dOL _Qo are nonzero. Differentiating Eq. (|4|), 
we obtain \dH,2/d0\g^o, = (|a|/2)[if,2(0°)V-ffc2(90^, 
where all parameters can be extracted experimentally. 
The temperature dependence of \dHc2/dO\Q^Qo is plot- 
ted in the inset of Fig. for all three superlattices. 
For a thin film with strong Pauli paramagnetic effect, 
the temperature dependence of \dHc2/d6\Q^Qo for the 
n = 5 superlattice can be expected, i.e. a small but 
increasing slope as the temperature is lowered, proba- 
bly dominated by an increasing Hp with decreasing tem- 
perature. However, for n = 4 and n = 3 superlattices, 
in particular n = 3, the behaviour is strikingly differ- 



ent in that \dHc2/dO\Q^Qo is a lot larger near Tc and 
shows an opposite trend in the temperature evolution. 
Since Hp{T) should increase upon cooling, the suppres- 
sion of \dHc2/dO\0^Qo at low temperatures can only be 
attributed to a rapid decrease of k. Our results therefore 
provide evidence that an FFLO-like inhomogeneous state 
is approached at low temperatures. 

We now provide a mechanism to explain the enhance- 
ment of \dHc2/dO\0^Qo as the thickness of the film is re- 
duced. In view of the fact that diverges near Tc, this 
exceptional sensitivity of \dHc2/d6\Q^Qo to layer thick- 
ness d is surprising, since d <C is satisfied for all three 
superlattices studied near Tc. Instead, it appears that 
the orbital pair breaking is a lot more effective in the 
n = 3 superlattice. To explain our data, we point out 
the importance of space inversion symmetry. In a sys- 
tem where space inversion symmetry is broken, it has 
been shown that the Pauli effect can be substantially re- 
duced for all magnetic field orientations [21] . Here, the 
inversion symmetry is locally broken at the top and the 
bottom 'interface' CeCoIns layers at the immediate prox- 
imity to YbCoIns block layers. Although the inner layers 
might also break the local inversion symmetry, the effect 
has been theoretically shown to be much weaker [TB]. 
Thus, with the reduction of n, the fraction of the non- 
centrosymmetric interface layers increases rapidly (Fig. 
[5|. We argue that local inversion symmetry breaking is 
responsible for the enhancement of \dHc2/d6\Q^Qo with 
decreasing n through the weakening of the Pauli effect. 

In the presence of the local inversion symmetry break- 
ing, coupled with the fact that cerium has a large atomic 
number, the Rashba-type spin-orbit coupling is expected 
to be strong [3], giving rise to FS branches where the de- 
gree of splitting is proportional to the coupling strength 
aR (Fig. [5|. In this case, because of the violation of spin 
and momentum conservations along with the large effec- 
tive mass, the coupling between a non-centrosymmetric 
layer with Rashba-split FS and a centrosymmetic layer 
where the FS is intact is expected to be weak, leading 
to an overall reduction of the interlayer hopping integral 
tc. In the limit where aR ^ tc^ interesting physics such 
as the enhancement of spin susceptibility have been pre- 
dicted by Ref. [16] . Therefore, by combining the prospect 
of manipulating both the time reversal and space inver- 
sion symmetries, these superlattices offer a new avenue 
for detailed study of exotic superconducting phases, e.g. 
a helical superconducting phase and novel vortex 
physics [23 . 

In summary, we have measured and analysed Hc2 (0) of 
CeCoIn5(n)/YbCoIn5(5) superlattices for n = 3,4 and 5. 
We have shown that the rounding of the cusp at ^ = 0° 
at low temperatures can be interpreted as a signature 
that the system approaches an FFLO-like phase. Fur- 
thermore, the drastic difference in the temperature de- 
pendence of Hc2{0) for different n can be understood by 
considering the importance of local inversion symmetry 
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breaking, the effect of which is much stronger in the n = 3 
superlattice. These artificial heavy fermion superlattices 
where the effects of Pauh paramagnetism and Rashba 
interaction entangled thus offer a new playground for ex- 
ploring exotic superconducting phases. 
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